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SYMBOLS 


c 

local speed of sound 

c p 

pressure coefficient 

E 

total energy 

H 

enthalpy 

MACH 

Mach number 

Mt 

maximum tip Mach number due to the rotation of rotor 

V 

pressure 

q 

absolute velocity in the inertial frame of reference 

Q r 

relative velocity in the rotor-fixed frame of reference 

s 

entropy 

SMACH 

sectional Mach number 

u 

forward flight velocity of the rotor 

x/c 

chordwise relative position on a profile 

1 

ratio of specific heats 

V 

advance ratio 

Q 

angular velocity of the rotor 

d 

partial differential operator 


blade azimuthal position 

P 

density 

e c 

collective pitch angle 
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SUMMARY 


A new method is presented for calculating the quasi-steady transonic flow over a lifting or 
nonlifting rotor blade in both hover and forward flight by using Euler equations. The approach 
is to solve the Euler equations in a rotor-fixed frame of reference using a finite volume method. 
A computer program was developed and was then verified by comparison with wind-tunnel 
data. In all cases considered, good agreement was found with available experimental data. 

INTRODUCTION 

The computational prediction of aerodynamic forces on three-dimensional lifting rotor 
blades is the first step toward the computerized design of advanced rotors. It is apparent that 
computational methods will play a major role in the design and analysis of future rotary-wing 
vehicles (ref. 1). Several prediction methods (ranging from transonic small disturbance (refs. 2- 
3), through transonic full-potential (refs. 4-10), to Euler (refs. 11-12) techniques) for calculat- 
ing transonic flow over a rotor blade have recently been proposed. However, the development 
of more accurate, efficient and robust methods will certainly continue until the viscous flow 
field around an entire rotorcraft can be economically calculated. 

In this paper, a new solution method is described for the numerical solution of the Euler 
equations as applied to helicopter rotors. The method solves the three-dimensional Euler equa- 
tions in a rotor-fixed frame of reference using a finite volume method. A computer program 
based on this method is designated as TFAR3 in the Transonic Flow Analysis for Rotors (TFAR) 
series of advanced computer codes developed at NASA Ames Research Center (refs. 5-6). The 
TFAR3 code can be used to calculate the flow past a lifting rotor, in both hover and forward 
flight, with the blade tip in a transonic flow regime. The geometry of the rotor blade can be 
quite general. The usual helicopter design parameters such as shaft angle, collective pitch an- 
gle, local twist angle, flapping angle, and variable tip shape are taken into account. 

The TFAR3 code had been verified by comparison with wind-tunnel data (refs. 13-14). 
To establish the usefulness of the code, computational results are presented for the following 
cases in the hierarchy of flow-field complexity. 

1. Nonlifting rotors in forward flight: computations were made for rotors with different 
blade tip shapes. 

2. Lifting rotors in hover: several comparisons were made for a two-bladed model rotor 
in hover at different tip speeds. 

3. Lifting rotors in forward flight: computations were made for a three-bladed lifting 
model rotor in forward flight. 

In all the cases considered, good agreement was found with published experimental data. 

This is Part 3 of a series of planned publications under the same general title "Tran- 
sonic Flow Analysis for Rotors." The author would like to thank Dr. Frank Caradonna and 


Dr. Chee Tung of the U. S. Army Aviation Research and Technology Activity, Aeroflightdy- 
namics Directorate, for their critical and fruitful discussions and their generous assistance with 
the rotor wake modeling. 


FLOW FIELD ANALYSIS PROBLEM 

The rotor flow field is acknowledged to contain very complicated three-dimensional, un- 
steady, and viscous flow phenomena. During each revolution, the rotor blade may enter a 
transonic flow regime on the advancing side and may face boundary-layer separation on the re- 
treating side. Moreover, the blade may encounter wakes shed from other blades with a complex 
blade-vortex interaction. Because of these complexities, the entire problem cannot be modeled 
with current computer capabilities. In this paper, we will focus on the advancing-blade side of 
the rotor flow field, particularly on the transonic effects. 

While the potential flow solutions have proved extremely useful for transonic flows with 
shock waves of moderate strength, they are limited by an inherent inability to correctly predict 
stronger shock jumps. Also, the roll-up of the near wake and the convection of tip vortices are 
not easily modeled because a complex wake tracking is required. In this section we describe 
the flow field analysis problem in terms of Euler equations and discuss the boundary conditions 
at the outer computational boundaries. 

Euler Equations in a Rotor-Fixed Reference Frame 

Let the rotor-fixed frame of reference move with a constant linear velocity U and a constant 
angular velocity Q with respect to an inertial frame of reference. The relative velocity V of a 
point fixed in the moving frame of reference is 

V= U+n xr 


where r is the position vector. Let the ambient density p^, speed of sound c a c , and the reference 
chord length be unity. The nondimensional Euler equations in this rotor-fixed frame of reference 
are found to be 

Pt + V(p9 r ) = 0 

(pq r )t + V(pg r g r + p) + 2pQ x q r + pQ x V + pV t = 0 

(pE) t + V(pHq r ) + p(q r +V)V t = 0 ( 1) 

where q T is the velocity in the moving frame, p the density, p the pressure, E the total energy, 
and H the total enthalpy. The absolute velocity in the inertial frame of reference is given by 
q = q r + V. 


For a perfect gas, 


and 


E ( 7 - l)p + 2 Qt ' Qr 2 V ' V 


H = E+ - 

P 


2 



where 7 is the ratio of the specific heats and is equal to 1.4 for air. 

The tangential flow boundary condition on the rotor blade surface is 

q r ■ n = 0 

where n is the unit vector normal to the surface. At the outer computational boundaries, non- 
reflection boundary conditions as detailed in the next subsection are imposed. 

Far-Field Boundary Condition 

The boundary conditions at the far-field computational boundaries are derived by intro- 
ducing Riemann invariants for one-dimensional flow normal to the boundaries. 

The characteristic form of the one-dimensional Euler equations can be written as 
pc[dt+ (<?r • n±c)dn](q • n) ±[d t + (q r ■ n± c)d n ]p + pc(Q x q ■ n) = 0 

[d t + {q r ■ n)d n ]S = 0 

where c is the speed of sound, n is the outward unit normal vector to the boundary, and S is the 
entropy. 

For local isentropic flow, let a thermodynamic quantity a(p) be defined by 



then the first equation, after dividing through by pc, is reduced to be 

[3t + (q r • n± c)d n )(q ■ n ± a) + Q x q ■ n = 0 


The quantities 


i? ± = q ■ n ± a + 


Q. x q ■ n 
q T ■ n ± c 


( r ■ n) 


which are called the Riemann invariants, are constant for the forward and backward character- 
istics. The invariants are propagated with speeds q r • n±c along the characteristics respectively. 
Let subscripts 00 and e denote the free-stream values and the values extrapolated from the in- 
terior adjacent to the boundary. For 7-law gas, considering subsonic flow, the free-stream and 
extrapolated Riemann invariants Roo and R e at the outflow computational boundary are found 


to be 


Roc ~ 9oo ■ n 


1 ^ + / 
7 - 1 


and 


2c e 

R e = q e - n+ + 9 

Of - 1 
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where f and g are given by 


/ = 


fl x goo - n 

Qroo ' 71 Cqo 


(r ■ n) 


and 


9 = 


£2 x q e • n 
Qre • n+ c e 


(r • n) 


correspondingly, and r is the position vector. These may be combined to give 


q-n= j(Re + Roc- f -g) 
and 

(l — 1 ) 

C= - L ^- L (R c -Roo + f-g) 

The quantities q ■ n and c are actually specified at the computational boundary. The other 
components of velocity are then extrapolated from the interior through the following formula 


q = q e + (q ■ n- q e ■ ri)n 


The entropy S is also extrapolated from the interior S = S e for the outflow boundary. The 
density, energy and pressure can then be calculated from the entropy and the speed of sound. 
This procedure avoids specifying the pressure at the downstream boundary where it does not 
necessarily recover its free-stream value because of the trailing wake. 

The conditions at the inflow, the inner radial and the outer radial computational boundaries 
can be derived in a similar manner. 

For a computational region large enough to enclose the entire flow field about the whole 
rotor and its wake, the above boundary value problem is complete. For a computational region 
covering only a portion of the rotor disk, the flow effect caused by the other rotor blades and 
the flow outside the finite computational region have to be incorporated into the boundary con- 
ditions in some way. There are several methods of doing this. In the present work, this effect 
will simply be supplied by the inflow angle correction technique (refs. 7,15). 

The quasi-steady problem is defined as the solution of the flow equations (1) with all the 
time derivative terms set equal to zero. 

NUMERICAL SOLUTION PROCEDURE 

Recently, efficient methods (refs. 16-19) for solving the Euler equations for steady flows 
have been developed for fixed-wing applications. Because the aim of this study is to develop 
a method for design work, the concern is for the easy applicability of the method to complex 
configurations. One possible approach to meet this need is the finite volume method. The 
Euler solver reported in references 19 and 20 has been adapted and modified for rotary-wing 
applications. 
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Finite Volume Scheme 


The integral form of the Euler equation (1) can be written as 

a ‘IIL wdT+ IL FnU+ IIL GdT=0 <2) 

Here, W is the conserved quantity, F the flux term, G the forcing term, dr the volume element, 
ds the surface element, and n the outward unit normal vector. 

The discretization procedure follows the method of lines in decoupling the approximation 
of the spacial and temporal terms. First, the computational domain in physical space is divided 
into hexahedral cells. Equation (2) is then approximated in each cell separately. Let the depen- 
dent variables be described at each center point (x, j , k) of the cells. A coupled set of ordinary 
differential equations is then obtained 

4-idrW) + Q+ T = 0 (3) 

at 

where Q is the flux term and T is the forcing term. The cell volume dr can be determined as 
the sum of six tetrahedrons. Each of them can be evaluated by a simple formula as given in 
reference 21. 

The flux term is given by 


QiJ,k 'y ^ 
celLsides 

where Aij,k denotes the area of the cell face between the points (i, j , A;) and (i + 1 , j , k). The 
face area can be determined as one-half of the vector product of the cell face diagonals. The 
value of Fij,k at the cell face is taken as the average of F at the points (t, j, k ) and (i + 1 , 
;, k). The use of centered differences ensures the scheme of second order accuracy in spatial 
discretization, provided that the mesh is smooth enough. 

At the blade surface, the only contribution to the flux balance comes from the pressure. 
The normal pressure gradient at the wall can be estimated in terms of quantities which can be 
determined from the interior solution as given in reference 22. The pressure on the blade is 
then estimated by extrapolation from the pressure at the adjacent cell centers. 

Dissipative Terms 

To inhibit odd/even point decoupling and to prevent the appearance of wiggles in the re- 
gions where pressure gradient is large, it is necessary to add artificial dissipative terms. Namely, 
equation (3) is replaced by 

(4) 


— (drW) + Q + T- D = 0 
dt 
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where D is the dissipation term. An effective form for D is a blend of second and fourth order 
differences with coefficients depending upon the second order difference of the local pressure. 
Specifically, let v be defined by 

_ I Pi+ijjt - 2 Pij,k + | 

~ : 

Pt'+l,;',* + 2 PiJ,k + Pi-\J,k 

This quantity is second order except in regions containing a large pressure gradient. The term 
D<j,k is constructed so that it is thud order in smooth regions of flow and first order in the 
vicinity of the shock wave. Dijjk has the form 

Dij,k = (D x + D y + D z )uij t k 

where u; is p for the continuity equation, pq for the momentum equations, pE for the energy 
equation, and D x > Dy and D z are central difference operators in i-, j- and /^-directions respec- 
tively. Specifically, D x is defined by 

D x Uij t k = d,+i/2j,jk — d-i-\f 2 j,k 


, _ d T i+\/2J,k f (2) a 

“«+l/2 J,k ^ '- € i+\/2J,k^ lU3 'J£ 


e «'+t /2 J,k^x w i-l J,*] 


and A x is the forward difference operator given by 


Ar WiJ,k ~ w i+ 1 J,k ~ W ij,k 


The coefficient 

€< il\/2j,k = « (2) ' max(v M j,k,Vij,k) 
is made proportional to v . The other coefficient is given by 

^i+l/2j,k = max( 0 .,k ( 4) - %\ /2Jtk ) 

The fourth order differences provide background dissipation throughout the computational 
domain but are switched off in the neighborhood of the shock waves. The typical values of the 
constants ac ( 2) and « (4) are 1/4 and 1/256. 

Runge-Kutta Scheme 

If the cell volume dr is constant in time, the equation (4) has the form 

W t + R(W) = W t + P(W) + D(W) =0 (5) 
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where the residual R(W) consists of convective part P(W ) and dissipative part D(W). A 
class of Runge-Kutta four-stage time-stepping scheme can be employed to solve this system of 
ordinary differential equations. It can be written as 

W {0) = W n 

W (1) = W (0) -aiAffl (0) 
w ( 2 ) = w ( 0 ) _ a2 Atfl(D 

( o) 

W {3) = w W _ a3 A tR {2) 

W (4) _ w ( 0 ) _ tR (3) 

= w (4) 

where At is the time step, W n and the value of W at the beginning and end of the n-th 
time step, and in the (/c+l)-th stage R (k) = P(W (k) ) + D( W (0) ) . The dissipation is evaluated 
once and then frozen after the first stage in each time step calculation. The choice of coefficients 
determines the characteristics of the scheme. For single grid calculations, a good choice of the 
coefficients is given by ai = .6, 0.2 = .6, <*3 = 1., and <24 = 1. 

Residual Smoothing 

One way to extend the stability range of the explicit scheme is to perform the implicit 
residual smoothing. This process extends the domain of dependence and thus increases the 
stability range of the scheme. 

The residual smoothing is applied in product form 

( 1 - rud 2 x ) • ( 1 - r yd]) • ( 1 - r\ k d 2 z )Rij,k = Rij,k 

where, Rij t k is the residual before smoothing and Rij,t is the new residual. The coefficients are 
given by 77 , = 1 , tj, = 1 , and Tjjt = 1 . Because it is only necessary to solve a sequence of tridiagonal 
equations for scalar variables, this scheme requires a relatively small amount of computational 
effort per time step. Comparatively, other implicit schemes usually need to solve a coupling 
system with a much more costly block-tridiagonal solver. It is well known that a Runge-Kutta 
scheme with properly tuned residual smoothing can be shown to be unconditionally stable based 
on linear stability analysis (ref. 20 ). However, in practice, the fastest convergence to a steady 
state solution is realized with a small time step. 

Local Time Stepping 

For a steady state problem, the convergence of the scheme can be accelerated by using 
local time step. Namely, the solution advances in time with a time step size dictated by the 
local stability limit. This allows faster signal propagation on coarse mesh regions and thus 
faster convergence of the overall scheme. For single grid calculations, the use of local time 
step of Courant-Friedrichs-Lewy (CFL) number 4.0 gives rise to a good convergence rate for 
most cases. 
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Wake Modeling 

For lifting rotor calculations, the wake effect and the blade motion are important. Because 
at present only a single blade of a rotor is contained in the finite domain of the computer code, 
a practical method is required to account for the far-field wake effect and other blade effects. A 
basic approach to solve this problem is to divide the entire flow field into two parts. The blade 
and its near wake will be modeled in the present code direcdy. The rest of the flow field which 
lies outside of the computational domain will be modeled through a set of inflow angles along 
the blade span at each azimuth position. This set of inflow angles can be furnished through 
auxiliary helicopter wake analysis codes. Two wake prediction codes are employed in the 
present study. For hover prediction, a linearized lifting surface code as reported in reference 23 
is used. For forward flight prediction, a lifting line code as described in reference 24 is utilized. 
This technique was originally reported in reference 15 and had been employed in several recent 
papers (refs. 6, 8, 10). 


RESULTS 

Since only the coordinates in physical space of each grid point of the computational region 
are needed in the finite volume scheme, the body aligned mesh can be generated by several 
methods. The present computer code for a single rotor blade uses a simple C-type grid, which 
can be generated by a combination of stretching and conformal mappings similar to that in 
reference 5. Good solutions may be obtained by employing a grid of 129 points in stream wise, 
33 points in normal, and 33 points in radial directions. The grid comprises 81 points on the 
airfoil surface and 21 points on the span of the rotor blade. A typical calculation needs about 
600 time steppings to obtain a good solution and will take less than 15 min of CPU time on the 
Cray-XMP computer at Ames Research Center. For the same case, a multigrid calculation will 
take about 100 cycles to get a comparably accurate solution and will take about 5 min of CPU 
time on the same Cray-XMP computer. 

The present calculations are mainly focused on transonic flows that occur on the 
advancing-blade side of a rotor in forward flight. All these studies have been performed with 
rigid blades in order to separate the aerodynamic and aeroelastic effects. A number of results 
are presented for the flow problems that are typically encountered in rotor flow calculations. 

Nonlifting Rotors in Forward Flight 

The first set of calculations is concerned with the fundamental case of nonlifting rotors 
for which there is no wake effect. Calculations were performed for model rotor blades with 
straight-tip and swept-back tip planforms. Figure 1 shows the planform of an ONERA nonlift- 
ing straight-tip model rotor blade (ref. 13). The surface pressure transducers are located at three 
spanwise stations of .855, .892, and .946. In figures 2a-d, the computed surface pressure distri- 
butions at these stations of the above blade are plotted and compared with ONERA measured 
data at an azimuth position i/> of 90° for moderate (0.4) to high (0.55) advance ratios /i. The 
tip Mach number Mr due to rotation for all the cases is 0.6. The results are quite satisfactory. 
Figure 3 shows the planform of the same ONERA model blade but with swept-back 30° blade 
tip. Figure 4 shows similar results for the model blade with this second blade tip at a tip Mach 
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number of 0.6288 and advance ratio of 0.5. The effect of this swept-back tip is clearly seen as 
a significant reduction of the suction peak over what occurs for a corresponding unswept case 
(fig. 2c). 


Lifting Rotors in Hover 

Calculations were made for a rigid two-bladed model rotor (ref. 14) in hover at various 
tip speeds. This model rotor employs an untapered NACA 0012 airfoil section and has an 
aspect ratio of 6. Calculations were performed for blades with a collective pitch angle 6 C of 8 0 
and at three different tip Mach numbers of .612, .794, and .877. Here, the wake effects were 
calculated through a lifting surface code for hover prediction (ref. 23). Figures 5a-c show the 
comparison of calculated and measured surface pressure distributions at three spanwise stations 
of .80, .89, and .96 for the above three tip speeds. The overall agreement of surface pressure 
distributions is excellent for these three different cases. The locations of shock waves and the 
lower surface pressure distributions are predicted very well. The predictions of the surface 
pressure distributions at the inner two spanwise stations, 0.5 and .68 (which are not shown in 
this paper), are not as good as those at the outer spanwise stations. The disagreement may 
result from the presence of two vortex lines in the near field which are not currently taken into 
account. 


Lifting Rotors in Forward Flight 

For forward lifting flight calculation, the wake effects were obtained through a lifting line 
code (ref. 24). Computed results have been compared with experimental data for an ONERA 
three-bladed lifting model rotor at various high-speed forward flight conditions (ref. 13). Fig- 
ure 6 shows the comparison of calculated and measured surface pressure distributions at three 
spanwise stations of .852, .902, and .951 near the blade tip. The azimuthal angle is 90° . The 
tip Mach number is 0.6288 and the advance ratio /i is 0.3872. The overall agreement of surface 
pressure distributions is good. The discrepancies seen on the lower surface are presently not 
understood. However, similar comparisons have been found with other codes (refs. 6, 10). 

CONCLUDING REMARKS 

A new solution method for calculating the inviscid, rotational, quasi-steady transonic flow 
about helicopter rotors has been constructed. This method solves the conservation form of the 
Euler equations in a rotor-fixed frame of reference by a finite volume method. Based on this 
method, a three-dimensional Euler code was developed and validated. It has been shown that 
the code can be employed to calculate the flow past a lifting rotor in hover as well as in forward 
flight with the blade tip in a transonic flow regime. The geometry of the rotor blade can be 
quite general. A number of cases have been solved by using this new computer code, and the 
results were found to be accurate. 
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The tip vortex formation might be modeled by using this Euler code, and this possibility 
will be studied in the future. Other future works include incorporating a known vorticity field 
into the Euler calculation, and extending the current computer code to solve the Navier-Stokes 
equations for viscous flow calculations. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, California 94035-1000, January 26, 1990 


10 



REFERENCES 


I Davis, S. S.; and Chang, I-Chung: The Critical Role of Computational Fluid Dynamics 
in Rotary-Wing Aerodynamics. AIAA Paper 86-0336, Jan. 1986. 

2 Caradonna, F. X.; and Isom, M. P.: Numerical Calculations of Unsteady Transonic Po- 
tential Flow over Helicopter Rotor Blades. AIAA Jnl., vol. 14, Apr. 1976, pp. 482-488. 

3 Chattot, J. J.: Calculation of Three-Dimensional Unsteady Transonic Flows Past Heli- 
copter Blades. NASA TP-1721, 1980. 

4 Arieli, R.; and Tauber, M. E.: Computation of Subsonic and Transonic Flow about Lifting 
Rotor Blades. AIAA Paper 79-1667, Aug. 1979. 

5 Chang, I-Chung: Transonic Flow Analysis for Rotors - Part I. Three-Dimensional, Quasi- 
Steady, Full-Potential Calculation. NASA TP-2375, 1984. 

6 Chang, I-Chung: Transonic Flow Analysis for Rotors - Part II. Three-Dimensional, Un- 
steady, Full-Potential Calculation. NASA TP-2375, 1985. 

7 Chang, I-Chung; and Tung, C.: Numerical Solution of the Full-Potential Equation for 
Rotors and Oblique Wings Using a New Wake Model. AIAA Paper 85-0268, Jan. 1985. 

8 Tung, C.; and Chang, I-Chung: Rotor Transonic Computation with Wake Effect. Fourth 
International Conference on Applied Numerical Modeling, Taiwan, Rep. of China, Dec. 1984. 

9 Sankar, N. L.; and Prichard, D.: Solution of Transonic Flow Past Rotor Blades Using the 
Conservative Full-Potential Equation. AIAA Paper 85-5012, Oct. 1985. 

10 Strawn, R. C.; and Caradonna, F. X.: Numerical Modeling of Rotor Flows With a Con- 
servative Form of the Full-Potential Equations. AIAA Paper 86-0079, Jan. 1986. 

II Roberts, T. W.; and Murman, E. M.: Solution Method for a Hovering Helicopter Rotor 
Using the Euler Equations. AIAA Paper 85-0436, Jan. 1985. 

12 Sankar, N. L.; Wake, B. E.; and Lekoudis, S. G.: Solution of the Unsteady Euler Equa- 
tions for Fixed and Rotor Wing Configurations. J. Aircraft, vol. 23, no. 4, Apr. 1986, 
pp. 283-289. 

13 Philippe, J. J.; and Chattot, J. J.: Experimental and Theoretical Studies on Helicopter 
Blade Tips at ONERA. ONERA TP 1980-96, Sep. 1980. 

14 Caradonna, F. X.; and Tung, C.: Experimental and Analytical Studies of a Model Heli- 
copter Rotor in Hover. NASA TM-8 1232, Sep. 1981. 

15 Caradonna, F. X.; Desopper, A.; and Tung, C.: Finite-Difference Modeling of Rotor 
Flow Including Wake Effects. 8th European Rotorcraft Forum, Aix-en-Provence, France, Aug. 
1982. 

16 Briley, W. R.; and MacDonald, H.: Solution of the Three Dimensional Compressible 
Navier-Stokes Equations by an Implicit Technique. Proc. 4th International Conference on 
Numerical Methods in Fluid Dynamics, 1974. 

17 Beam, R. M.; and Warming, R. F.: An Implicit Finite Difference Algorithm for Hy- 
perbolic Systems in Conservation-Law Form. J. Comp. Phys., vol 22, no. 1, Sep. 1976, 
pp. 87-110. 

18 MacCormack, R. W.: A Numerical Method for Solving the Equations of Compressible 
Viscous Flow. AIAA Paper 81-110, Jan. 1981. 


11 


19 Jameson, A.; Schmidt, W.; and Turkel, E.: Numerical Solution of the Euler Equations 
by Finite Volume Methods Using Runge-Kutta Time Stepping Schemes. AIAA Paper 81-1259 
Jun. 1981. 

20 Jameson, A.; and Baker, T. J.: Solution of the Euler Equations for Complex Configura- 
tions. AIAA Paper 83-1929, Jul. 1983. 

21 Kordulla, W.; and Vinokur, M.: Efficient Computation of Volume in Flow Predictions. 
AIAA Jnl., vol. 21, no. 6, Jun. 1983, pp. 917-918. 

^Rizzi, A.: Numerical Implementation of Solid Body Boundary Conditions for the Euler 
Equations. ZAMM, vol. 58, no. 7, Jul. 1978, pp.301-304. 

23 Summa, J.M.: Advanced Rotor Analysis Methods for the Aerodynamics of Vortex- 
Blade Interactions in Hover. Vertica, vol. 9, no. 4, 1985, pp. 331-343. 

24 Johnson, W.: Development of a Comprehensive Analysis for Rotorcraft 1 . Rotor Model 
and Wake Analysis. Vertica, vol. 5, no. 2, 1981, pp. 99-129. 


12 



r- <N CO 



Figure 1. - Planform geometry of an ONERA nonlifting model rotor blade with straight 
tip. 


Mj = 0.6, JU = 0.4, it = 90° 



SMACH = 0.755 0.777 0.809 

Figure 2. - Comparison of computed and measured surface pressure distributions for a 
nonlifting straight-tip blade in forward flight, (a) Advance ratio of .4. 
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M T = 0.6, ii = 0.55, \p = 90 
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Figure 2. - Concluded, (d) Advance ratio of .55. 
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Figure 3. - Planform geometry of an ONERA nonlifting model rotor blade with swept- 
back 30° tip. 
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M t = 0.6288, fi = 0.50, \j/ = 90° 



SMACH =0.852 0.884 0.913 

Figure 4. - Comparison of computed and measured surface pressure distributions for 
nonlifting swept-back-tip blade in forward flight. 
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M T = 0.877, 0 C = 8 
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STATION= 0.8000 



0.8900 



0.9600 


MACH =0.7016 
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Figure 5. - Concluded, (c) Tip Mach number of .877. 
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M T = 0.6288, n = 0.3872, »// = 90° 



X/C X/C X/C 


STATION= 0.8520 0.9020 0.9510 

MACH =0.7792 0.8107 0.8415 

Figure 6. - Comparison of computed and measured surface pressure distributions for a 
lifting three-bladed rotor in forward flight. 
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